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Abstract 

We study current fluctuations in lattice gases in the macroscopic limit extending the 
dynamic approach for density fluctuations developed in previous articles. More precisely, 
we establish a large deviation principle for a space-time fluctuation j of the empirical 
current with a rate functional We then estimate the probability of a fluctuation 

of the average current over a large time interval; this probability can be obtained by 
solving a variational problem for the functional X. We discuss several possible scenarios, 
interpreted as dynamical phase transitions, for this variational problem. They actually 
occur in specific models. We finally discuss the time reversal properties of I and derive a 
fluctuation relationship akin to the Gallavotti-Cohen theorem for the entropy production. 
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1. Introduction 

1.1. Thermodynamic junctionals for non equilibrium systems 

In equilibrium statistical mechanics there is a well defined relationship, established by 
Boltzmann, between the probability of a state and its entropy This fact was exploited 
by Einstein to study thermodynamic fluctuations. So far it does not exist a theory of 
irreversible processes of the same generality as equilibrium statistical mechanics and pre- 
sumably it cannot exist. While in equilibrium the Gibbs distribution provides all the 
information and no equation of motion has to be solved, the dynamics plays the major 
role in non equilibrium. 

When we are out of equilibrium, for example in a stationary state of a system in contact 
with two reservoirs, even if the system is in a local equilibrium state so that it is possible to 
define the local thermodynamic variables e.g. density or magnetization, it is not completely 
clear how to define the thermodynamic potentials like the entropy or the free energy. One 
possibility, adopting the Boltzmann-Einstein point of view, is to use fluctuation theory 
to define their non equilibrium analogs. In fact in this way extensive functionals can 
be obtained although not necessarily simply additive due to the presence of long range 
correlations which seem to be a rather generic feature of non equilibrium systems. This 
possibility has been pursued in recent years leading to a considerable number of interesting 
results. One can recognize two main lines. The first, as well exemplified by the work of 
Derrida, Lebowitz and Speer [12-14], consists in exact calculations in specific models of 
stochastic lattice gases. The second is based on a macroscopic dynamical approach for 
Markovian microscopic evolutions, of which stochastic lattice gases are a main example, 
that leads to some general variational principles. We introduced this approach in [3, 4] 
and developed it in [5,6]. Both approaches have been very effective and of course give the 
same results when a comparison is possible. 

Let us recall the Boltzmann-Einstein theory of equilibrium thermodynamic fluctua- 
tions, as described for example in [23]. The main principle is that the probability of a 
fluctuation in a macroscopic region of fixed volume V is 

P oc exp{VAS/k} (1.1) 

where AS is the variation of the specific entropy calculated along a reversible transforma- 
tion creating the fluctuation and k is the Boltzmann constant. Eq. (|1.1|) was derived by 
Einstein [16] simply by inverting the Boltzmann relationship between entropy and prob- 
ability. He considered (jl.l)) as a phenomenological definition of the probability of a state. 
Einstein theory refers to fluctuations for equilibrium states, that is for systems isolated or 
in contact with reservoirs characterized by the same chemical potentials. When in contact 
with reservoirs AS is the variation of the total entropy (system + reservoirs) which for 
fluctuations of constant volume and temperature is equal to —AF/T, that is minus the 
variation of the free energy of the system divided by the temperature. 

We consider a stationary nonequilibrium state (SNS), namely, due to external fields 
and/or different chemical potentials at the boundaries, there is a flow of physical quanti- 
ties, such as heat, electric charge, chemical substances, across the system. To start with 
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it is not always clear that a closed macroscopic dynamical description is possible. If the 
system can be described by a hydrodynamic equation, a fact which can be rigorously 
established in stochastic lattice gases, a reasonable goal is to find an explicit connection 
between the thermodynamic potentials and the dynamical macroscopic properties like 
transport coefficients. As we discussed in [3-6], the study of large fluctuations provides 
such a connection. It leads in fact to a dynamical theory of the free energy which is 
shown to satisfy a Hamilton- Jacobi equation in infinitely many variables requiring as in- 
put the transport coefficients. In the case of homogeneous equilibrium states the solution 
of the Hamilton- Jacobi equation is easily found, and the equilibrium free energy is recov- 
ered together with the well known fluctuation-dissipation relationship, widely used in the 
physical and physical-chemical literature. On the other hand in SNS the Hamilton- Jacobi 
equation is hard to solve. There are few one-dimensional models where it reduces to a non 
linear ordinary differential equation which, even if it cannot be solved explicitly, leads to 
the important conclusion that the non equilibrium free energy is a non local functional of 
the thermodynamic variables. This implies that correlations over macroscopic scales are 
present. The existence of long range correlations is probably a generic feature of SNS and 
more generally of situations where the dynamics is not invariant under time reversal [2]. 
As a consequence if we divide a system into subsystems the free energy is not necessarily 
simply additive. 

Besides the definition of thermodynamic potentials, in a dynamical setting a typical 
question one may ask is the following: what is the most probable trajectory followed by the 
system in the spontaneous emergence of a fluctuation or in its relaxation to an equilibrium 
or a stationary state? To answer this question one first derives a generalization of the 
Boltzmann-Einstein formula from which the most probable trajectory can be calculated 
by solving a variational principle. The free energy is then related to the logarithm of 
the probability of such a trajectory and satisfies the Hamilton- Jacobi equation associated 
to this variational principle. For equilibrium states and small fluctuations an answer to 
this type of questions was given by Onsager and Machlup in 1953 [27]. The Onsager- 
Machlup theory gives the following result under the assumption of time reversibility of 
the microscopic dynamics: the most probable creation and relaxation trajectories of a 
fluctuation are one the time reversal of the other. As we show in [3,4], for SNS the 
Onsager- Machlup relationship has to be modified in the following way: the spontaneous 
emergence of a macroscopic fluctuation takes place most likely following a trajectory which 
can be characterized in terms of the time reversed process. 

1.2. Macroscopic dynamics and large fluctuations 

We consider many-body systems in the limit of infinitely many degrees of freedom. 
Microscopically we assume that the evolution is described by a Markov process X T which 
represents the state of the system at time r. This hypothesis probably is not so re- 
strictive because also the dynamics of Hamiltonian systems interacting with thermostats 
finally is reduced to the analysis of a Markov process, see e.g. [15]. To be more precise 
X T represents the set of variables necessary to specify the state of the microscopic con- 
stituents interacting among themselves and with the reservoirs. The SNS is described by 
a stationary, i.e. invariant with respect to time shifts, probability distribution P st over 
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the trajectories of X T . We denote by p the invariant measure of the process X T . The 
measure p is a probability on the configuration space and for each fixed time r we have 
P st {X T = u) = 

We assume that the system admits a macroscopic description in terms of density fields 
which are the local thermodynamic variables pi. The usual macroscopic interpretation of 
Markovianity is that the time derivatives of the thermodynamic variables p\ at a given 
instant of time depend only on the p^s and the affinities (thermodynamic forces) |^ at the 
same instant, recall that F is the free energy As we discussed in [6], for non equilibrium 
systems, the affinities, defined as the derivative of the non equilibrium free energy, do 
not determine the macroscopic evolution of the variables p$. There is an additional non 
dissipative term which however does not modify the rate of approach to the stationary 
state. 

For simplicity of notation we assume that there is only one thermodynamic variable p 
e.g. the local density. For conservative systems the evolution of the field p = p(t, u), where 
t and u are the macroscopic time and space coordinates, is then given by the continuity 
equation 

d tP = V- [±D(p)Vp-x(p)E] = -V-J(p) (1.2) 

where D(p) is the diffusion matrix, x{p) t ne mobility and E the external field. Finally the 
interaction with the reservoirs appears as boundary conditions to be imposed on solutions 
of (jl.2j) . We shall denote by p = p{u) the unique stationary solution of (jl.2j) . i.e. p is the 
typical density profile for the SNS. 

This equation derives from the underlying microscopic dynamics through an appro- 
priate scaling limit in which the microscopic time and space coordinates r, x are rescaled 
diffusively: t = t/N 2 , u = x/N where N is the linear size of the system so that the 
number of degrees of freedom is proportional to N d . The hydrodynamic equation (jl.2j) 
represents a law of large numbers with respect to the probability measure P st conditioned 
on an initial state X Q . This conditional probability will be denoted by Px - The initial 
conditions for ()1.2|) are determined by X®. Of course many microscopic configurations 
give rise to the same value of p(0, u). In general p = p(t, u) is the limit of the local density 
ir N (X T ). 

The free energy F(p), defined as a functional of the density profile p = p(u), gives the 
asymptotic probability of density fluctuations for the invariant measure p. More precisely 

p(n N (X) « p) ~ exp { - N d F(p)} (1.3) 

where d is the dimensionality of the system, ttn(X) ~ p means closeness in some metric 
and ~ denotes logarithmic equivalence as iV — > oo. In the above formula we omitted the 
dependence on the temperature since it does not play any role in our analysis; we also 
normalized F so that F(p) = 0. 

In the same way, the behavior of space time fluctuations can be described as follows. 
The probability that the evolution of the random variable ttn(X t ) deviates from the 
solution of the hydrodynamic equation and is close to some trajectory p(t) is exponentially 
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small and of the form 



P st (n N {X N 2 t ) « p(t), t E M 2 ]) ~ exp { - iV d [F(p(^)) + ^lO 5 )] } (1.4) 



where T{p) is a functional which vanishes if p(t) is a solution of ()1.2|) and F(p(ti)) is the 
free energy cost to produce the initial density profile p(t\). Therefore T{p) represents the 
extra cost necessary to follow the trajectory p(t) in the time interval [ti, t?\. Equation (jl.4j) 
is a dynamical generalization of the Boltzmann-Einstein formula, we shall refer to it as 
the dynamical large deviation principle with dynamicalrate functional T . For stochastic 
lattice e shown in [4], the functional T can be calculated explicitly. 

To determine the most probable trajectory followed by the system in the spontaneous 
creation of a fluctuation, we consider the following physical situation. The system is 
macroscopically in the stationary state p at t — — oo but at t = we find it in the state 
p. According to (jl.4j) the most probable trajectory is the one that minimizes T among 
all trajectories p(t) connecting p to p in the time interval [— oo, 0]. As shown in [3,4] this 
minimization problem gives the non equilibrium free energy, i.e. 



To this variational principle it is naturally associated a Hamilton- Jacobi equation which 
plays a crucial role in the analysis developed in [3-6]. We emphasize that the func- 
tional J 7 , hence the corresponding Hamilton- Jacobi equation for F, is determined by the 
macroscopic transport coefficients D(p) and x{p)i which are experimentally accessible, 
see e.g. [1]. We can thus regard (jl.5|) as a far reaching generalization of the fluctuation- 
dissipation theorem since it allows to express a static quantity like the free energy in terms 
of the dynamical macroscopic features of the system. 

1.3. Current fluctuation and related thermodynamic junctionals 

Beside the density, a very important observable is the current flux [11,18,25,28,29]. 
This quantity gives informations that cannot be recovered from the density because from 
a density trajectory we can determine the current trajectory only up to a divergence free 
vector field. We emphasize that this is due to the loss of information in the passage from 
the microscopic level to the macroscopic one. 

In the previous paper [7] we have introduced a Boltzmann-Einstein type formula for 
current fluctuations. This formula shows that the asymptotic probability, as the number 
of degrees of freedom increases, of observing a current fluctuation j on a space-time 
domain [0, T] x A can be described by a rate functional X[ 0j T](j)- In the present paper we 
develop the approach introduced in [7] and illustrate some relevant applications. 

To discuss the current fluctuations, we introduce a vector-valued observable J N \{X a , 
< a < t}) of the trajectory X which measures the local net flow of particles. As for the 
density, for stochastic lattice gases, we shall be able to derive a dynamical large deviations 
principle for the current. Recall that Px stands for the probability P st conditioned on 
the initial state Xq. Given a vector field j : [0, T] x A — ► M. d , we have 



F(p) = inf F^^p) 



(1.5) 



p 



P Xo (J N (X) « j(t, «)) ~ exp { - N d X [0 ,T](j)} 



(1.6) 
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where the rate functional is 



Wi) = lj o dt([j - J(p)U{ P r l \3 - J{P)\) (1-7) 
in which we recall that 

J(p) = -±D(p)Vp + X (p)E . 

Moreover, p = p(t, u) is obtained by solving the continuity equation d t p + V ■ j = with 
the initial condition p(0) = po associated to X . The rate functional vanishes if j = J(p), 
so that p solves (jl.2j) . This is the law of large numbers for the observable J N . Note 
that equation () 1.7(1 can be interpreted, in analogy to the classical Ohm's law, as the total 
energy dissipated in the time interval [0,T] by the extra current j — J{p)- 

The functional I describes the fluctuation properties of the current, the density and 
all observables related to them, as proved in Section El Among the many problems 
we can discuss within this theory, we study the fluctuations of the time average of the 
current J N over a large time interval. This is the question addressed in [9] in one space 
dimension by postulating an "additivity principle" which relates the fluctuation of the 
time averaged current in the whole system to the fluctuations in subsystems. We show 
that the probability of observing a given divergence free time average fluctuation J can 
be described by a functional $(J) which we characterize, in any dimension, in terms of a 
variational problem for the functional I[o,t] 

$(J)= lim infij [0iT] (j), (1.8) 

where the infimum is carried over all paths j = j{t,u) having time average J. The 
static additivity principle postulated in [9] gives the correct answer only under additional 
hypotheses which are not always satisfied. Let us denote by U the functional obtained by 
restricting the infimum in (jl.BJl to divergence free current paths j, i.e. 

U(J) = inf \{[J - J{p)l X {pY 1 [J - J{p)\) (1.9) 

where the infimum is carried out over all the density profiles p = p(u) satisfying the 
appropriate boundary conditions. From (jl.8)l and (jl.9j) it follows that $ < U. In one 
space dimension the functional U is the one introduced in [9]. 

There are cases in which $ = U and in Subsection 16. II below we give sufficient condi- 
tions on the transport coefficients D, x for the coincidence of $ and U. On the other hand, 
while $ is always convex the functional U may be non convex. In such a case U (J) under- 
estimates the probability of the fluctuation J. In [7] we interpreted the lack of convexity 
of U, and more generally the strict inequality $ < U, as a dynamical phase transition. In 
the present paper we investigate in more detail the occurence of this phenomenon. Let 
us denote by U** the convex envelope of U; then $ < U** and in Subsection 16.31 we give 
an example where U** < U . 
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We shall also consider the fluctuation of the time averaged current with periodic 
boundary conditions. In Subsection 16.21 we discuss the behavior of U and <3> under ap- 
propriate conditions on the transport coefficient and the external field. In particular 
we show that for the Kipnis-Marchioro-Presutti (KMP) model [20], which is defined by 
a harmonic chain with random exchange of energy between neighboring oscillators, we 
have U(J) = (l/2)J 2 /x(m) = (l/2)J 2 /m 2 , where m is the (conserved) total energy. In 
addition we show, for J large enough, $(J) < U(J). This inequality is obtained by 
constructing a suitable travelling wave current path whose cost is less than U(J). We 
mention that, by using the space-time approach introduced in [7], the possibility of tak- 
ing advantage of travelling waves has been first envisaged by Bodineau and Derrida [10] 
for the periodic simple exclusion process with external field. We refer to the discussion in 
Subsection 16.21 for a comparison between KMP and simple exclusion models. 

We study also the behavior of X and $ under time reversal and derive a fluctuation 
relationship akin to the Gallavotti-Cohen theorem for the entropy production [17,22,24]. 
In fact, we prove, in the present context of lattice gases, that the anti-symmetric part of 
$ is equal to the power produced by the external field and the reservoirs independently 
of the details of the model. From this relationship we derive a macroscopic version the 
fluctuation theorem for the entropy production. 

2. Microscopic model 

As the basic microscopic model we consider a stochastic lattice gas with a weak exter- 
nal field and particle reservoirs at the boundary. The process can be informally described 
as follows. We consider particles evolving on a finite domain. At each site, independently 
from the others, particles wait exponential times at the end of which one of them jumps 
to a neighboring site. Superimposed to this dynamics, at the boundary particles are 
created and annihilated at exponential times. More precisely, let A C R d be a smooth 
domain and set An = NA n 1* d . We consider a Markov process on the state space X An , 
where X is a subset of N, e.g. X = {0, 1} when an exclusion principle is imposed. The 
number of particles at site x G is denoted by r\ x G X and the whole configuration by 
r) G X An . The dynamics is specified by a continuous time Markov process on the state 
space X An with infinitesimal generator Ln = N 2 [L ^n + £&,iv] defined as follows: for 
functions / : X An -> R, 

Lo, N f(v) = \ E c *MW x ' V r))-m], 

\x-y\=l 

Wfo) = \ E { c *M[f(v x ' v v)-f(v)}+c y Av)[f(v v > x v)-f(v)}}- 

xeAjv,»/0Ajv 
\x-y\=l 

(2.1) 

Here \x\ stands for the usual Euclidean norm and a x,y r], x,y G Ajv, for the configuration 
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obtained from 77 by moving a particle from x to y: 

{Vz if z ^ x, y 

+ 1 if z = V 
rj x — l if z = x . 

If x £ A^r, y £" A^r, then <7 y,x ?7 is obtained from rj by creating a particle at x, while a 2- '' 5 '?] is 
obtained by annihilating a particle at x. Therefore the generator Lo,tv describes the bulk 
dynamics which preserves the total number of particles whereas L^n models the particle 
reservoirs at the boundary of An- Note that we already speeded up the microscopic time 
by N 2 in the definition of Ln, which corresponds to the diffusive scaling. 

Assume that the bulk rates c XtV , x, y £ A^v, satisfy the local detailed balance [24] with 
respect to a Gibbs measure defined by a Hamiltonian TC and in presence of an external 
vector field E = (Ei, . . . , Ed) smooth on the macroscopic scale. Likewise, assume that the 
boundary rates c XtV , c VtX , x £ Ajv, y £" A N , satisfy the local detailed balance with respect 
to 7i and in presence of a chemical potential Xo(y/N) smooth on the macroscopic scale. 

The above requirements are met by the following formal definitions. Fix a smooth 
function A : A — > R and a Hamiltonian H. Consider jump rates c° x y satisfying the 
detailed balance with respect to the Gibbs measure associated to 7i with free boundary 
conditions if x, y £ An, while if x £ An, y ^ An we add the chemical potential Xo{y/N): 

cl y (rj) = exp-{H(<T x «ri)-n(v)}cl x (cT*>y V ), x,yeA N ; 

cl y ( v ) = ex P -{H(a x >y V ) - H(rj) + Hy/N)}^^) , x £ A N , y <? A N . 

Note that we included the inverse temperature in the Hamiltonian Ti. Of course if r\ x = 
then c° xy (r)) = 0. 

Fix a smooth vector field E = (E\, . . . , Ed) ■ A — > R d and let 

c x , x+e M ■= e N ~ lE ^ cl x+ei (r,) , c x+ei ,M := e- N - lE ^ c° x+eu M , (2-2) 

where {ei, . . . , e^} stands for the canonical basis in R d . Namely, for large, by expanding 
the exponential, particles at site x feel a drift N~ 1 E(x/N). 

Typically, for a non equilibrium model, we would consider A as the rf-dimensional cube 
of side one, the system under a constant field E/N and a chemical potential A satisfying 
X (y/N) = 70 if the first coordinate of y is 0, X (y/N) = 71 if the first coordinate of y is 
N, imposing periodic boundary conditions in the other directions of A. 

By setting c XiV = if both x and y do not belong to An, we can rewrite the full 
generator L N as follows 

N 2 

L N f{ri) = — Yl ^Mlf{o^v)-m] (2-3) 

x,y& d 
\x—y\=\ 

We consider an initial condition 77 £ X An . The trajectory of the Markov process r)(t), 
t > 0, is an element on the path space -D(R + ; X An ), which consists of piecewise constant 
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paths with values in X Ajv . We shall denote by the probability measure on D (M+; X Ajv ) 
corresponding to the distribution of the process rj(t), t > with initial condition rj. 

Examples of stochastic lattices gases are the simple exclusion processes in which X = 
{0, 1}, Ti — and c° y (r/) = 7} x [l — %] and zero range processes in which X — N, Tl(r)) = 
J2a>J2i<k<ri(x)^ &9(k)i ^ or some function g : N — »■ K+ such that g(0) = 0, and c° xy (rf) = 
9(Vx)- 



3. Macroscopic description of lattice gases 

The empirical density tt n can be naturally defined as follows. To each microscopic 
configuration rj G X An we associate a macroscopic profile n N (u), u G A, by requiring that 
for any smooth function G : A — > R 



7T . 



G) = /d« 7r"(u) G(u) = ^ E G ( x / N )^ (3- 1 ) 



so that ir N (u) is the local density at the macroscopic point u = x/N in A. Of course 
n N (u) is really a sum of point masses at the points x/N with weight r) x /N d ; in the limit 
N — > oo it will however weakly converge to a "true" function 

The definition of the empirical current is slightly more complicated. Indeed it is not 
a function of the configuration 77 G X An but of the trajectory {^(t)}t>o £ -D(M+; X Ajv ). 
Given an oriented bond (x, y), let N x,y (t) be the number of particles that jumped from x 
to y in the time interval [0, £]. Here we adopt the convention that Af x,y (t) is the number 
of particles created at y due to the reservoir at x if x ^ A N , y G A N and that Af x,y (t) 
represents the number of particles that left the system at x by jumping to y if x G A^r, 
y ^ Ajy. The difference Q x,y (t) = J\f x,y (t) — M y,x {t) is the net number of particles flown 
across the bond {x,y} in the time interval [0, t]. Given a trajectory r)(s), < s < t, the 
instantaneous current across {x,y} is defined as dQ x,y /dt. This is a sum of 5-functions 
localized at the jump times with weight +1, resp. —1, if a particle jumped from x to y, 
resp. from y to x. 

For a given realization of the process r)(t) in £)(R + ; X Ajv ), we define the corre- 
sponding empirical current J N as follows. Let T > and pick a smooth vector field 
G = (Gi, . . . , G d ) defined on [0, T] x A. We then set 



{{J N ,G)) T = [ dt [ duG(t,u)-J N (t,u) 
Jo Ja 



= ]v4r E E / T ^ ^' X+ei W . (3-2) 

where • stands for the inner product in M d and we sum over all x such that either x G Ajy 
or x + ti G Atv. The empirical current J N is therefore a signed measure on ([0, T] x A) , 
while we recall that the empirical density is a positive measure on A. The normalization 
AM^+i) in flE) has been chosen so that the empirical current has a finite limit as A r — > 00. 
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The local conservation of the number of particles is expressed by 

Vx{t)-Vx(o)+ Yl Q x ' y (t) = o. 

y:\x-y\=l 

It gives the following continuity equation for the empirical density and current. Let G 
be a smooth function on a neighborhood of the closure of A. Denote by VatG the vector 
field whose coordinates are (V ^G) .(u) = N[G(u + Ci/N) — G(u)]. Then 

(ir N (T), G) - (n N (0), G) = ((J N , V N G)) T - f ± £ G(y/N)dQ^(t) 

Jo xeA N ,y&A N 
\x-y\=l 

The above equation can be formally stated as the continuity equation 

d t n N + V N -J N = (3.3) 

In particular, given the initial condition, the trajectory of the process, described by the 
empirical density n N can be completely recovered from the empirical current J N . 

We briefly discuss at the heuristic level the law of large numbers, as iV — > oo, for 
the empirical density and the empirical current. Details are given in Appendix [X] Fix 
a sequence of configurations 7] N and assume that its associated empirical measure ir N 
converges to p (u)du for some density profile po : A — > M + . Let us denote by p = p(t,u), 
J = J(t,u), the limiting values of ir (t, u), <J N (t,u), respectively. Here n (t,u) is the 
empirical density associated to the configuration rj (t) and <J N (t,u) has been defined in 

(J321). 

The microscopic relation ()3.3|) implies the continuity equation 

+ V • J = (3.4) 

To derive a closed evolution for p and J, we need to express the current J in terms of the 
density p. To simplify the exposition, we assume the process to be gradient: there exist 
local functions (rj), i = l,...,d, depending on the configuration r/ around 0, so that 
for any i = 1, • • • , d 

c° x , x+ei (v) ~ c° x+ei , x (v) = hg>{ V ) - h { X e M 

where h x is the function evaluated on the configuration r] translated by x. 

Denote by p>\ the infinite volume grand canonical ensemble relative to the Hamiltonian 
TC with chemical potential A. Choose the chemical potential A = A(p) so that fix[r] ] = p 
and define 

d ( Hp) = »m[ h f] . X {1) (P) ■= (1/2K( P ) K, ei + eg, >0 ] (3.5) 
We show in Appendix El that the current J can be expressed in terms of the density p as 

J=-~D(p)Vp + X (p)E=: J(p) (3.6) 
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where D and x are d x d diagonal matrices with entries Du(p) = j-d^'(p) and Xu(p) — 

x (l) (p)- 

For non gradient systems the diffusion matrix D and the mobility x are not in general 
diagonal. In such a situation D is given by a Green-Kubo formula [30, II. 2. 2] and x 
can be obtained by linear response theory [30, II.2.5]. These coefficients are related by 
Einstein relation D = R~ 1 x, where R is the compressibility: = Fq, in which F$ is 
the equilibrium free energy associated to the Hamiltonian 7i, [30]. 

To conclude the description of the evolution, it remains to examine the evolution at 
the boundary of A. We claim that the density is fixed there because we speeded up 
diffusively the non-conservative Glauber dynamics at the boundary: 

X(p(t,u)) = X (u) uedA (3.7) 

The macroscopic evolution of the density and the current is thus described by the 
equation 

" dtp + V • J = , ue A 
J=-\D{p)Vp + x {p)E, ueA 
X(p(t,u)) = Ao(u) , uedA 
p(0,-) = Po(-)- 

The stationary density profile p = p(u), u e A, is the stationary solution of the 
hydrodynamic equation, that is 

V • J(p(u)) = , u E A 
A(p(u)) = AoO) uedA, 

If we let the macroscopic time diverge, t — > oo, p(t) — > p and J(p(t)) converges to J(p), 
which is the current maintained by the stationary state. 

We next discuss the large deviation properties of the empirical current. More details 
are given in Appendix EI As before we consider a sequence of initial configuration rj 
such that the empirical density 7r N (r] N ) converges to some density profile p Q . We fix a 
smooth vector field j : [0, T] x A — >• M. d . The large deviation principle for the current 
states that 

F% N (j N (t,u)*j(t,u), (t,u) e [0,T] xA)~ exp{ - N d I [0>T] (j)} (3.8) 
where the rate functional X is 



1 



T 



%t](j) = - / dt{w) - j{p{t)%x{pit)r i w) - j( P (t))]) (3.9) 



in which p(t) = p(t, u) is obtained by solving the continuity equation 

d t p(t : u) + V -j(t } u) = 
p(0,u) = p (u) 



(3.10) 



10 



and J(p) is given by (|3.6I) . 

Of course there are compatibility conditions to be satisfied, for instance if we have 
chosen a j such that p(t, u) becomes negative for some (t, u) G [0, T] x A then X[ ,r](j) = 
+00. Notice that, even if not indicated explicitly in the notation, the rate functional X 
depends on the initial density profile po, through equation ()3.1()j) . 

We note that in the large deviation functional ()3.9|) the fluctuation of the density p(t) 
is determined by the current j(t). The large deviations properties of the density, which 
we described in [3-5] for non equilibrium stochastic lattice gases, can thus be deduced 
from the ones of the current, see Appendix El for the details. This is due to the fact that 
the continuity equation, as already remarked, holds exactly at the microscopic level, see 
(J3.3JI . On the other hand the constitutive equation (|3.6|) holds only in the limit A" — > 00 
when fluctuations can be neglected. 



4. Large deviation of the time averaged current 

We want to study the fluctuations of the time average of the empirical current over a 
large time interval [0, T]; the corresponding probability can be obtained from the space 
time large deviation principle ()3.8j) . Fix T > and a divergence free vector field J = J(u). 
We introduce the set of possible paths j of the current with time average J 

A T ,j = {j = j(t,u) : i J dtj{t,u) = J(u)} 

The condition of vanishing divergence on J is required by the local conservation of the 
number of particles. By the large deviations principle ()3.8|) , for T and A" large we have 



N 



£dt J N (t) « j) ~ exp { - N d T$(J)} (4.1) 



where the logarithmic equivalence is understood by sending first N — * 00 and then T — > 
00. In Subsection 16.41 below we shall show that for the zero range process the limits can 
be taken in the opposite order; we expect this to be true in general. The functional $ is 
given by 

®(J) = I™ inf izj&ntf) = inf .mf ^X[o,T](j) (4.2) 

By a standard sub-additivity argument we show that the limit T — > oo exists and coincides 
with the infimum in T. Indeed, given ji G At,j and ji £ As,j, we have 

Z[o,t+s] U) = %t] (ji) + Z[o,s\ (h) (4.3) 

where j is obtained by gluing j 1 and j 2 . Here we used the invariance of X under time 
shift and that j\ G At,j implies that the corresponding density p±, obtained by solving 
the continuity equation ()3.10jl . satisfies pi(0) = Pi(T). From the previous equation we 
get the sub-additivity property: 

. pf X[ ,r+s](j) < . inf X [0iT ](j) + . mf X [0>s] (j) 
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Even if the rate functional X depends on the initial density profile po, by taking the limit 
in (|4.2|) it is easy to show $ does not. 

We now prove that $ is a convex functional. Let < p < 1 and J = pJ\ + (1 — p) J 2 , 
we want to show that $(</) < p&( J\) + (1 — p)Q(J2)- By ([4.2jl . given e > we can find 
T > 0, ji e Apr,Ji, and j 2 e A(i_ p )t,j 2 so that 

> ^^O'O-e 

$ ( J 2) > 7^ I ^; ^[o,(i-p)t](J2)-£ 

By the same arguments used in ([4.3)1 . the path obtained by gluing ji with j 2 , denoted by 
j, is in the set At,j- Therefore, 

HJ) < ^[0,T](j) < P $(Jl) + (l-p)$(J 2 )+£ 

which proves the convexity of $. These arguments are standard in proving the existence 
and convexity of thermodynamic functions in statistical mechanics. 

We next study the variational problem on the right hand side of (14.2)1 . We begin by 
deriving an upper bound. Given p = p{u) and J = J(u), V ■ J = 0, let us introduce the 
functionals 

U{p,J) = \(J-J(P),X(P)- 1 [J-J(P)]) (4-4) 

U(J) = inf U(p, J) (4.5) 

p 

where the minimum in ()4.5)) is carried over all profiles p satisfying the boundary condition 
([3.711 and J(p) is given by ()3.6[1 . When J is constant, that is, in the one-dimensional case, 
the functional U is the one introduced in [9]. 
We claim that 

HJ) < U(J) . (4.6) 

The strategy to prove this bound is quite simple, see also [9]. Let p = p(J) be the density 
profile which minimizes the variational problem ([4.5)1 . Given the initial density profile p$, 
we choose some fixed time r > and a current j = j(u) which moves the density from 
Po to p in a time lag r, namely such that rV ■ j = po — p. We now construct the path 
j = j(t, u), (t, u) e [0, T] x A as follows 

{j if < t < r 
-^J if T <t<T-T 
-j if T-r <t<T 

The corresponding density p{t) is obtained by solving the continuity equation ()3.10p . i.e. 

( P0 + HP-P0) if 0<t<r 
p(t) = < p if r <t<T -t 

{ Po + ^(p-po) if T-r <t <T 
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It is straightforward to verify that j G At,j, as well as lirriT^oo -f 2T[o,T](j) = U(J). 

By the convexity of $( J) we can improve the upper bound ()4.6|) for free. Let us denote 
by U** the convex envelope of U, i.e. the largest convex functional below U. By taking 
the convex envelope in (|4.fij) we get 

$(J) < U**(J) (4.7) 

We next discuss a lower bound for the variational problem ([4.2)1 . We denote by W 
and [/ the same functionals as in (|4.4jl - (|4.5p . but now defined on the space of all currents 
without the conditions of vanishing divergence. Let also U** be the convex envelope of 
U. Let j G At, j- By the convexity of U** in the set of all currents, we get 

^Wl(j') = ^£dtU(p(t),j(t)) > ±£dtU(j(t)) 



1 



T 



> - y o dtu**(j(t)) > u**(j) 

which implies 

$(J) > £T*(J) (4.8) 

The upper and lower bounds (|4.7|) and ()4.8|) are different in general. For a divergence free 
J we have £/(</) = U(J) but since the convex envelopes are considered in different spaces, 
we only have U**(J) < U**(J). 

The derivation of the upper bound shows that our result differs from the one in [9] if 
U is not convex. Moreover, if $(</) < U(J), the optimal density path p in the variational 
problem (|4.2jl must be time dependent. 

We now examine how different behaviors of the solution to the variational problem 
()4.2|) reflect different dynamical regimes that we interpret as dynamical phase transitions. 
It is convenient to work in the time interval [— T, T] instead of [0,T]. We consider the 
system in the ensemble defined by conditioning on the event (2T)~ 1 J^ T dt J N it) = J with 
N and T large. The parameter J plays therefore the role of an intensive thermodynamic 
variable and the convexity of $ expresses a stability property with respect to variations 
of J. 

If $(</) = U(J) and the minimum for ()4.5|) is attained for p = p(J) we have a state 
analogous to a unique phase: by observing the system at any fixed time t we see, with 
probability converging to one as N, T — > oo, the density 7r N {t) ~ p(J) and the current 
J N (t) ~ J. 

When $(J) = U**(J) < U(J), we have a state analogous to a phase coexistence. 
Suppose for example J = pj x + (1 - p) J 2 and U(J) > U**(J) = pU(J x ) + (1 - p)U(J 2 ) 
for some p, J\, J 2 - The values p, J±, J 2 are determined by J and U. The density profile 
is then not determined, but rather we observe with probability p the profile p{J\) and 
with probability 1 — p the profile p{J 2 )- Actually there is a memory of initial condition: 
if we take p(—T) = p(T) = p(Ji) we will see the density p{J\) in the time intervals 
[-T,-(l-p)T] and [(1 — p)T, T], p(J 2 ) in [-(1 - p)T, (1 - p)T). 
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Consider now the case in which a minimizer for (|4.2|) is a function j(t) not constant in 
t. This is possible (an example will be given in Subsection 16. 2j) only when $(J) < U**(J). 
Suppose first that j(t) is periodic with period r and denote by p(t) the corresponding 
density Of course we have r _1 J^dt j(t) = J. In such a case we have in fact a one 
parameter family of minimizers which are obtained by a time shift a £ [0, r]. By choosing 
2T an integral multiple of r and p(—T) = p(a) for some a £ [0, r] then the empirical 
density in the conditional ensemble will follow the path p(t + a + T). This behavior is 
analogous to a non translation invariant state in equilibrium statistical mechanics, like 
a crystal. Finally if j(t) is time dependent and not periodic the corresponding state is 
analogous to a quasi-crystal. 

The asymptotic (j4.1j) can be formulated in terms of the Laplace transform of the 
empirical current as follows. For each divergence free, time independent, vector field 
A = A (it) we have 

where $*(A) is the Legendre transform of $(</): 

$*(A)= sup {(A, J) -$(,/)}, 

where the supremum is carried over all the divergence free vector fields J. It follows from 
(@~nj) that [/*<$*. 

We conclude this section deriving a variational expression for U*. Recall the definitions 
0131), (g3J) of 

Z7*(A) = su P {(A,J)-i([J- JWI.xCp)- 1 ^-^)])} 

= sup { - i([J - J(p) - *(P)A], x(p) - V " - X(P)A]> 
J,p 1 2 

+ i(A,x(p)A) + (A,J(p))} 

To compute the supremum over J we decompose the vector field J(p) + x(p)A as follows 
J(p) + X(P)A = x(p) + [J(p) + x(p) (A - VV)] (4.10) 

where ^ solves 

V ■ (x(p)V^) = V ■ ( J(p) + x(p)A) u £ A 
-0(w) = u E dA 

Since the second term in the decomposition ()4.10j) is divergence free we get 

U*(\) = sup {-I(V^x(pW) + ^(A,x(p)A} + (A, J(p)}} (4.11) 
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where the supremum is over all density profiles p satisfying Fq(p(u)) = \o(u), u G dA. 

5. Time— reversal and Gallavotti— Cohen symmetry 

In this Section we discuss the properties of the rate functional for the current under 
time reversal. We also show that the functional $, which measures the probability of 
deviations of the time averaged current, satisfies a fluctuation theorem analogous to the 
Gallavotti-Cohen symmetry. 

5.1. Time-reversal properties of the rate functional 

In the previous Sections we discussed a large deviation principle given a fixed initial 
condition r] N associated to a density profile p , i.e. n N (r] N ) — > p . Now we consider 
instead the stationary process, namely the initial condition is distributed according to 
the invariant measure p, which is defined by ^ p N ' {rj)L^f{rj) = for any observable 
/ : X An — > R; recall the generator L]y has been defined in (12.3)1 . As discussed in the 
Introduction, the large deviations of the empirical density under the distribution p N are 
described by the non equilibrium free energy F, i.e., 

p N (7r N ^p) ~exp{-iV d F(p)} (5.1) 

In [3,4] we show that the functional F, which for equilibrium states is trivially related 
to the free energy, can be characterized by a variational problem on the dynamical rate 
functional for the density T introduced in (jl.4|) . see also (|A.15J) . To this variational 
problem is associated a Hamilton- Jacobi equation which plays a crucial role. 

In order to analyze the large deviations properties of the stationary process, since the 
initial condition is not fixed, it is natural to consider the joint fluctuations of the empirical 
density and current. We have 

^ N {n N ^p,J N ^j te[-T,T})~ex V {-N d g hTiT] (p,j)} (5.2) 

Here P^v, a probability measure on the space D(WL; X Ajv ), is the stationary process. 
Of course, fluctuations of the density and of the current are not independent since the 
continuity equation d t p + V • j = must be satisfied. Therefore the large deviation 
functional is 

( F(p(-T))+I hT , T] (j) if d t p + V- j = 
^ T] (P,J) = ^ +QO otherwige (5.3) 

where X has been introduced in (|3.9J) . If we are interested only in the current fluctuations 
in the stationary process we get the appropriate rate functional by projecting ([5.3)1 . 

mf{F(p )+l h T,T](j)} ■ 

Po 

Let us denote by the adjoint of the generator Ln 1)2.3)1 with respect to the invariant 
measure p N . We call the process generated by L%, which is still Markovian, the adjoint 
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process. We remark that the invariant measure of the adjoint process is again /i . Given 
a path r] G -D(M;X Ajv ) its time reversed is naturally defined as [1^77] (t) = r)(—t). The 
stationary adjoint process, that we denote by P^T, is the time reversal of P^jv, i-e. we 
have P^ a = P^v o We extend the definition of the time reversal operator $ to the 
current as (t) = —j(—t). Note that the current j changes sign under time-reversal. 
Then 

(n N * P ,J N *jte [-T, T]) = P^T (tt" « *p, J N « ,y te [-T, T]) 
At the level of large deviations this implies 

G[-T,T\(j>,j)=Gt- TtTl (0p,0j) (5.4) 

where G?L TT i ls the large deviation functional for the adjoint process. 

The relationship (|5.4JI has far reaching consequences. We next show that it implies a 
fluctuation dissipation relation for the current. We assume that the adjoint process has 
a dynamical large deviations principle of the same form as (|5.H|) with 1 replaced by X a 
where 

ii-T, T] (j) = \ fjt m) - n P {t))U{p{t)Y l \j{t) - J a (p(tm , 

in which J a {p) is the typical value of the current of the adjoint process. We divide both 
sides of (|5.4|) by 2T and take the limit T — > 0. By using (|5.3|) we get 

(^ dtp ) = \( J ~ W'Xip)-^ - J(p)\) - \{j + J a {plx{p)- X [j + J a (p)]) 

recalling that d t p + V ■ j = 0, this is equivalent to 

-(^, v • j) = -{ J{p) + rip), X {p)- X j) + \(J{p) + J a (p), xipYVip) - J a (p)}) 

which has to be satisfied for any p and j. By using that SF/5p vanishes at the boundary 
of A, see [4], we can integrate by parts the left hand side above and get 

J{p) + J a (p) = ~x(p)V 5 -f (5.5) 

dp 

(J(p), X{PVJ{P)) = (J a (p), x(p)"V a (p)> (5.6) 

Equation (|5.5J) is a fluctuation dissipation for the current analogous to the one for the 
density discussed in [4]. It also extends the relationships between currents and thermo- 
dynamic forces, see e.g. [26], to a non equilibrium setting. By plugging (|5.5|) into (|5.6j) we 
also get another derivation of the Hamilton- Jacobi equation mentioned before, i.e. 
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Let us now consider the variational problem (|4.2J1 as well as the same problem for the 
functional X a , we denote by $ a the corresponding functional. From (|5.4|) we get 

$(J) = $°(-J) (5.7) 

For reversible process this symmetry states that the functional $ is even. 

Let us consider a path j(t), t £ [— T, T] such that (2T) _1 J_ T dt j(t) = J for some 
divergence free vector field J. Recalling (jH.fi j) and that D(p)x(p)~ 1 = Fq(p) we have 

x{pY 1 J{p) = -\vFi{p) + E 

Since Fq(p(u)) = \q(u), u £ <9A, by developing the square in (j3.9j) and integrating by 
parts we get 

^G[-T, T] (P,J) = ^G[-t,t\{#pM ~ 2<J, E) + J dE A J ■ n (5.8) 

where dH is the surface measure on <9A and n is the outward normal to A. In particular 
this relation implies that if p,j is an optimal path for the variational problem defining 
$( J) then $p, dj is an optimal path for the variational problem defining $(— J). 
By taking the limit T — ► oo in (j5.8j) we get 

$(J) - $(-J) = $(J) - $ a (J) = -2{J,E) + / rfSAoJ-n (5.9) 

which is a Gallavotti-Cohen type symmetry in our space time dependent setup for macro- 
scopic observables. Note that the right hand side of (j5.9j) is the power produced by the 
external field and the boundary reservoirs (recall E is the external field and Ao the chem- 
ical potential of the boundary reservoirs). 

5.2. Entropy production 

Recall that we denote by ¥^ N the stationary state and by P^ a its time reversed, i.e. 
the stationary adjoint process. In the context of Markov processes the Gallavotti-Cohen 
observable is defined as 

1 dP^' a 
^ (T) = -2TiV^ l0g ^Pf 

where the subscript means that we consider both distributions in the time interval [-T, T]. 
We introduced the factor IT N d in order to discuss the asymptotic N, T — > oo. As 
discussed in [24, §2.4], Wn(T) can be interpreted as the microscopic production of the 
Gibbs entropy. For N fixed and T — > oo the functional Wn satisfies a large deviation 
principle with rate function /jv namely, 

¥^ N (W N (T)^q) ~ exp{ -2TN d f N (q)} (5.11) 



(5.10) 

-T.Tl 
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In [17,22,24] it is shown that Wn satisfies the Gallavotti-Cohen symmetry, which states 
that the odd part of Jn is linear with a universal coefficient: /at(?) — /at( — q) = —Q- 

An elementary computation, analogous to the one in [24], shows that, for the stochastic 
lattice gases as introduced in Section |2l we can express the functional Wn in terms of the 
empirical current. More precisely, we have 

j=l x x&Ajy 

y^A N 

where the summation is carried over all x such that either x G or x + ej G A^v- The 
previous equation can be understood as an entropy balance. Indeed, in the right hand 
side the first term is, for N large, the difference of the non equilibrium free energy at 
times T and — T, the second is the difference of the energy and the third is the work done 
by the external field E and the boundary reservoirs. Therefore, Wn can be interpreted 
as the total entropy produced by the system in the time interval [— T, T] . 

Recalling the definition of the empirical current J N , we can rewrite the above equation 

as 

«wn = h{-U l °*^wk +n(n{T)) - nH - T)) 

+ 2((j»,E}) hT>T] - J2 My/mQ'-'(\-T,T])\ (5 ' 12) 

We emphasize that, while the empirical current is a vector in M d , W N (t) is a scalar. 
From the previous expression it follows that, for any 5 > 

lim lim P^v ( W N (T) - 2(E, J{p)) - f dE X J(p) ■ h > s) = (5.13) 

T-»ooiV— oo V J dA J 

where we recall that J(p) = — (l/2)D(p) Vp + x(p)E is the typical current. 

We note that as T — * oo we can neglected the first line on the r.h.s. of (J5.12|) because 
it is a boundary term. We thus define 

Wn(T) = ^{2((j\E)) hT , T] ~^- d J2 Uy/N)Q^[-T,T]} (5.14) 

xgAjv 

y£Ajv 

which satisfies, as T —> oo with N fixed, the large deviation estimate (|5.11|) with the same 
rate function f N - 

On the other hand, since Wn(T) is a function of the empirical current, we can apply 
the large deviation principle (|5.2|) . We then get, by taking first the limit N — > oo and 
then T — > oo, 

F^(W N (T) k?) ~ exp { - 2TN d f(q)} (5.15) 
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where the rate function / can be expressed in terms of the functional X namely 



f{q) = lim inf —^ { _ TT[ {j) 



in which we introduced the set of currents 

fT 

2T 



2 \ dt{j(t),E) - I dt f dEX j(t)-n 

. J-T J-T JdA 



q\ (5.16) 



where we recall <i£ is the surface measure on <9A and n is the outward normal to A. 

Finally, since E and Ao are time independent we can take the time average of the 
empirical current in ()5.16|) . Recalling ()4.2|) . we get 



f(q) = inf $(J) (5.17) 



where 



\j:2(J,E)- X J-n = q, V- J = o) (5.18) 

where we inserted the condition V- J = because other things do not happen. The content 
of the variational problem ()5.17j) is to look for, among all possible currents, the best one 
to have a fixed entropy production. It is straightforward to verify that the symmetry 
f!5.9|) implies the classical Gallavotti-Cohen symmetry for the limiting functional /, i.e. 
f(q) — f(—q) = —q- On the other hand, if d > 1, equation (|5.9j) is more general than the 
classical Gallavotti-Cohen symmetry. 

In the above argument we first took the limit N — > oo and next T — » oo, but we 
expect that these limits could be taken in any order. In particular these would imply 
limjv^oo /aK?) = /(<?)• m Section IB~H we prove that this is the case for the zero range 
process. We finally note that in the one-dimensional case, setting A = [0, 1], we can easily 
solve (JSHTJ). We get 

f{q) = K2(£)-[A (l)-Ao(0)] 



6. Dynamical phase transitions: examples 

As we have discussed in Section 0J we always have the following inequalities 

U**(J) < $(J) < U**(J) < U(J) (6.1) 

for any divergence free J. A natural question is when the above inequalities are strict 
and when are equalities, in particular when $ = U. As discussed in Section HJ the 
strict inequality $(</) < U(J) is a dynamical phase transition on the ensemble defined 
by conditioning on the event in which the time average current equals J. In this Section 
we discuss several examples which show that different scenarios actually do take place in 
concrete models. As we have shown in Section El the macroscopic behavior (including the 
probability of large fluctuations) of the system is determined by the transport coefficients 
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D(p) and x(p)- 111 this Section we consider these as given functions and discuss the 
properties of the variational problem defining $. Specific choices of D and x correspond 
to well studied microscopic models, such as the simple excision processes, the zero range 
process, and the KMP model [20]. 

In Subsection 16.11 we find sufficient conditions on D and \ implying $ = U. In 
Subsection 16.21 we discuss periodic boundary conditions: under appropriate conditions 
on the transport coefficient and the external field we show that the minimizer for the 
variational problem defining U is obtained when p is constant (in space). Moreover, by 
considering travelling waves, we find for J large a better (space-time dependent) strategy 
so that $ < U. These conditions hold in particular for the KMP model with no external 
field. Moreover, for the exclusion process with sufficiently large external field, we show 
that there exists a travelling wave path of current whose cost is strictly less than the 
constant (in time and space) one. This was first observed in [10]. In Subsection 16.31 we 
give an example where U is non convex which implies $ < U. Finally, in Subsection 16.41 
we compute the Legendre transform of U for the one dimensional zero range process in the 
presence of external field. As a by product, we show that the macroscopic limit N j oo 
and T j oo can be interchanged. 

6.1. A sufficient condition for $ = U 

We consider the case when the matrices D(p) and x(p) are multiple of the identity, i.e., 
there are strictly positive scalar functions still denoted by D(p), x{p)i so that D{p) i: j = 
D(p)5ij, x(p)i,j = x{p)^i,ji h3 = 1) • • - id. We denote derivatives with a superscript. Let 
us first consider the case with no external field, i.e. E = 0, we shall prove that if 

D(p)x"(p) < D\p)x\p) for any p (6.2) 

then $ = U. In this case U is necessarily convex. 
Moreover we show that if 

D(p) X "(p) = D'(p)x'(p) for any p (6.3) 

then we have $ = U for any external field E. We mention that under the condition ()6.3j) . 
as shown in [8, §7], also the non equilibrium free energy F can be computed explicitly 
and it is a local functional. 

Condition (|6.2|) is satisfied for the symmetric simple exclusion process, where D = 1 
and x(p) = p(l — p)i P £ [0, 1]. Condition (j6.3J) is satisfied for the zero range model, 
where D(p) = ^f'(p) and x(p) = ^{p) fo r some strictly increasing function \I/ : IR + — > WL + . 
Condition ()6.3|) is also satisfied for the non interacting Ginzburg-Landau model, where, 
p G R, D(p) is an arbitrary strictly positive function and x(p) is constant. 

Let us consider first the case when E — and condition (|fj.2jl holds. In view of (j3.9|) 
and (|4.4|) . to prove that $ = U it is enough to show that for each j = j(t, u) G At,j, i.e., 
such that T" 1 J Q dt j(t) = J, and p(t) such that d t p(t) + V • j(t) = we have 

^Z [0 ,T](j) = ^£dtU(p(t),j(t)) > U(J) (6.4) 
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Instead of p we introduce a new variable a so that a = d(p) := f p dp'D(p'). Condition 
(J6.2j) is then equivalent to the concavity of the function X(a) := x(d -1 (a)) where d~ x is 
the inverse function of d. We introduce the functional 



V(a,j) := U{d-\a),j) = ~(j + ~Va, ^-Jj + ^Va]) 

where we used f|3.6p . We claim that the functional V is jointly convex in (a, j). Let us 
first show that this implies the lower bound ()6.4|) . We have 

±£dtU(p(t)J(t)) = ±£dtV(a{t),j{t)) 

> V(i^ dta{t),^J dtj{tj) > miV(a,J) = m£U{p,J) = U{J) 

in which we used the convexity of V in the second step and j e *4t,j in the third. 
To prove that V is jointly convex we write 

V(a,j) = sup V a (a, j) , V a (a,j) := (j + Jv«,a) - ^<a,X(a)a) 

and the supremum is taken over all smooth vector fields a = a(u) on A. Since X(a) is 
concave, for each fixed a the functional V a is jointly convex. The claim follows. 

In the case with non vanishing E, we can use the same argument, but the functional 
V is given by 

V(a,j) := U(d~\a),j) = ±(j + l -Va - X{a)E,^-[j + Uja - X(a)E] ) 

Condition (|fj.H|) is equivalent to X"(a) = 0; in this case we can easily show, as before, 
that V is jointly convex. 

6.2. Periodic boundary conditions 

In this subsection we consider the case when A = T, the one dimensional torus of side 
length one and constant external field E. If there is no external field, E = 0, then it is 
an equilibrium model; non equilibrium if E ^ 0. 

In this case we have the possibility of constructing a space time path (p(t,u), j(t,u)) 
of density and current in the form of a travelling wave for the variational problem (14.21) 
defining the functional $. We shall see that, under some assumptions on the transport 
coefficient D(p),x(p), this strategy for large J is more convenient than taking a density 
path p constant in time, so that $(J) < U(J). 

In the context of periodic boundary conditions, the proof that $ = U presented in the 
previous section applies if D is constant. In other words, we have that 3>(J) = U(J) for 
all J provided D is constant, \ is concave and E = 0. 
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Let M. m (Y) the convex set of positive functions p on the torus T such that J^du p(u) 
m; we call m the mass of p. In this context, 



U{J) = inf - / du 
p 2 Jo 



o x(pH) 

where the infimum is carried over Ai m (T) and J(p) was defined in (|3.6[) . For each v G R, 
let : R -> R + be defined by 

t /n ( 1 {J + ^[p(n) -m] - J(p)} 2 

^ J = mf- / du^ ; , , 6.5 

where the infimum is carried over Ai m (T). 
We claim that 

$ < V v (6.6) 

for each v. Indeed, consider a profile po in A^ m (T). Let T = w" 1 and set p(t,u) = 
p (u — vt), j(t,u) = J + t'[Po(^ — tv) — m ] i n the time interval [0, T]. An elementary 
computation shows that the continuity equation holds and that the time average over the 
time interval [0, T] of j(-,u) is equal J. In particular, 



*(J) < ^J*dtU{p{t),3{t)). 



On the other hand, it is easy to show by periodicity that the right hand side is equal to 

1 f 1 du {J + v [pa(u) - m] - J(p )} 2 



x(Po(u)) 

Optimizing over the profile po, we conclude the proof of (J6.6j) . 

Fix a mass m, an external field E and a current J. If J 2 /x + E 2 X is a convex function 
then 

and the optimal profile for the variational problem defining U(J) is the constant profile 
p(u) = m. In particular if 1/ X and E 2 X are convex functions then (|6.7|) holds for any J 
so that C/ is trivially convex. 

Indeed, fix a mass m, a current J and an external field E. For any profile p in A4 m (T), 

= r^v^m + (V |(i/2)ww|2 







X(p) Jo X(p) Jo 



because the cross term vanishes upon integration Here d(p) is the function introduced in 
the previous Subsection. We have 

( l du U- E X(P)} 2 = f l du \jP E? x {p)] - 2EJ > + E 2 X (m) - 2EJ 

Jo X{p) Jo L x(p) J " X(m) 
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where we used Jensen inequality, the convexity of J 2 /x + F 2 x, and p G M. m (Y). 
Therefore, for all profiles p in JA m (T), 



1 d \J ~ J{p)f > {J-E X (m)Y 



x(p) xM 

Since the cost of the constant profile p(u) = m is (1/2){J — £'x(m)} 2 /x(m), ()6.7|) is 
proven. 

For the KMP model [8,20] we have D(p) = 1 and x(p) = P 2 - Since x an d 1/x are 
convex functions, the assumptions for ()6.7|) are satisfied for any external field .E. For the 
simple exclusion process we have D(p) = 1 and x{p) = p(l — p)- m t ne case °f no external 
field, = 0, since x is concave and x _1 is convex, it satisfies both the hypotheses for 
$ = U and the ones for (|6.7|) : hence $(J) = (1/2) J 2 /m(l — m). If 7^ the assumptions 
for (|6.7|) holds only if \E/J\ is small enough. 

Fix a mass m,eGK and take the external field E = eJ. Assume that 

(l-e 2 X (m) 2 )x"(m) > (6.8) 

We claim that there exists w in R such that 

y w j(J) {l-e X (m)} 2 

h ^— K 2 X M • (6 - 9) 

where ^ has been defined in (|6.5j) . 

Fix a mass m, a current J, an external field E = eJ and take u = wJ. For .M m (T), 
by expanding the square we get that 

1 du {J + WJ[P ~ m] + ^ W ^ ~ Ex ^ (6 10) 

x(p) 

P I dn {l + Mr-rn)-ex(p)y 1 /"^.^Ip)] 2 



x(p) 4 7 x(p) 

because the cross term vanishes. Expand the square on the first integral. Let F(r) = 
F w ,m{ r ) be the smooth function defined by 

F(r) = i ^ ^- - 2e + e 2 X (r) • 

An elementary computation shows that 

F"(m) = —^-^^2x{m) 2 w 2 - 4x{m)x'{m)w + 2x'(m) 2 - x(m)x"(m) + eV'(m)x(m) 3 j . 

Let to = x'{ m )/x( m )- F° r this choice F"(m) < 0. In particular, we can choose a non 
constant profile p(u) in _A/f. m (T) close to m such that F"(p(u)) < for every u. Hence, by 
Jensen inequality, the coefficient of J 2 in (J6.10j) is strictly less than 

{1 - ex(m)} 2 _ 
X(m) 
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The statement follows. 

Fix e, assume that l/x + e 2 x is convex and ()6.8|) holds. Take E = eJ. Then, by (|6.6)1 . 
(16. 7|) and ()6.9jl we have < U(J) for all sufficiently large currents J. The KMP model 
satifies the above requirements if |e| < 1/m 2 . 

The exclusion process with external field E satisfies (16.8)) for \E\ > \ J\/m(l — m) but, 
for these values of E and J the function 1/x + e 2 x is not convex so we do not know 
whether U is given by ()6.7|) . By ()6.9|) we have that, for large current J, there exists a 
travelling wave whose cost is strictly less than the one of the constant profile p(u) = m. 
This is not enough to prove the strict inequality < U(J). 

We conclude by giving, for the KMP process with no external field, a description of 
the mechanism for the strict inequality $ < U in terms of the power necessary, according 
to ([A.lOp . to substain a time average current J. To get U(J) = (1/2) J 2 /m 2 we switch on 
a constant external field equal to J/m 2 which provides exactly the power U(J). On the 
other hand we can impose a time average current J by imposing a space time dependent 
external field of the type F(u — vt); the corresponding density and current paths are 
then travelling waves. By exploiting the convexity of p 2 we have shown that the second 
strategy, at least for J large, requires less power. 

6. 3. An example with non convex U 

We discuss here a special choice of the macroscopic transport coefficients D and x f° r 
which the functional U defined in (14.5)1 is not convex. In particular the upper bound (j4.7j) 
with the convex envelope differs from ([4.6)1 . 

We take d — 1, A = (0,1), E — 0, D(p) = 1, and x(p) a smooth function with 
x(0) = = (accordingly the density satisfies < p < 1) such that there exist 

0<i<B<l, 1st for which x{p) = e ~ ip if A < p < B. Furthermore we take the 
equilibrium boundary conditions p(0) = p(l) = p = (A + B)/2. We show that, for a 
suitable choice of the parameters A, B,£, there are Ji < J2 so that U"(J) < for any 

Je{J 1 ,J 2 ). 

Although we did not construct explicitly a microscopic lattice gas model in which the 
macroscopic transport coefficients meet the above requirements, we believe it would be 
possible to exhibit a model which has the same qualitative behavior. 

For d = 1, D = 1, E = 0, and J el, the Euler-Lagrange equation for the variational 
problem (j4.5j) defining the functional U(J) is 

2 P(n) = -^H)[ J "4 P(n) ]' (6.11) 
p(0)=Po, p(l)=Pi 

For the above choice of x an d of the boundary conditions, provided 4\J\ < B — A, a 
solution of (|6.11)) is given by 

~ / \ 2, coshLMu - 1/2)1 

as can be easily verified. Note indeed that A < pj < B since we assumed 4| J\ < B — A. 
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We shall later prove that, under the above conditions, (j6.12j) is the unique solution of 
the boundary value problem (|6.11|) . A simple computation then gives 



u(J) = ^y tanh y 

Let F(z) := ztanhz and z* be the unique positive root of z~ x = tanhz. Then F"(z) = 
2(1 -tanh 2 z)(l -ztanhz) < for z > z*. Hence U"(J) < if |J| G (2z*/£, (£ - A)/4). 
This interval is not empty provided I is chosen large enough. 

To show that (|6.12j) is the unique solution of the boundary value problem let us 

first prove that, given J ^ 0, any solution of (jb.ll)) satisfies the a priori bound \p'\ < 2\J\. 
Since p(0) = p(l), we can exclude the possibility that > 2\J\ for every u G [0, 1]. By 

the continuity of p', it is therefore enough to prove that |p'( w )l 2| J\ for every u G [0, 1]. 
Suppose conversely that there exists u* G [0, 1] such that p'{u*) = 2 J. Then, by the 
uniqueness of the Cauchy problem \p" = ^y[J 2 ~ iP' 2 ]> P( u *) = P*> P'( u *) — 2 J, we 
would get that the solution of (j6.11j) is p{u) = p(u*) + 2J{u — u*). Since this function 
does not satisfy the boundary conditions in (JbMl)) we find the desired contradiction. 

Since 4| J\ < B — A, the a priori bound \p'\ < 2\J\ implies that any solution p of (jfi.llj) 
satisfies A < p < B . For such values we have that x'{p)/x{p) — Uniqueness of the 
solution to (jSHH) can then be easily proven by explicit computations. 

6.4- Zero range processes 

In this section we consider the so-called one-dimensional zero-range processes which 
models a non-linear diffusion of lattice gases [19] under constant external field E. The 
model is described by positive integer- valued variables r] x representing the number of par- 
ticles at site x. The particles jump with rates (1/2)^(77^)) exp{E/N} to right, (l/2)g(rj x )) 
exp{— E/N} to the left, respectively. The function g(k) is such that g(k + 1) — g(k) > a 
for some a > and g(0) = 0. We assume that our system interacts with particle reservoirs 
at the sites and iV whose activity is given by ipo, cp±. 

The generator of this Markov process is given by ()2.3|) with Ajv = {1, • • • , N} and 



^x,x+l 



(v) = gin.) e E/N , c x+1 , x (v) = g(v* + i) e~ E/N 



for 1 < x < N — 1. Moreover, at the boundary, 

c,M= Voe E/N , c lfi ( V )=g( Vl )e- E / N , 

c N , N+ i{r]) = g(r} N ) e E/N , c N+ltN (rj) = ip% e~ E/N . 

Let V(u) = Eu and let <Pn(x) be the solution of 

for 1 < x < N and with boundary condition <Pn(Q) = (po, <£n(N + 1) = <fii- Here A N 
stands for the discrete Laplacian. 
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The invariant measure pn is the grand-canonical measure pn = ELeA obtained 
by taking the product of the marginal distributions 

( i \ 1 VN(x) k 



Z{ip N {x)) g(l) ■ ■ ■ g(k) 

where Z(ip) = Y^k>of k /idW ' ' '9{fy] * s ^ ne normalization constant. Let R(<p) = ipZ'(ip)/ 
Z(ip) and denote by \P its inverse function. For this process, the hydrodynamic equation 
(11.21) and the large deviations principle (jl.4j) can be obtained with D(p) = ^f'(p), x{p) = 
see [4]. Moreover, Fq(p) = log^(p). 
Let first show how, for this model, it is possible to solve explicitly the variational 
problem (J4.ilj) for the Laplace transform of the total current. We note that the case E — 
has already been solved in [9] . Since we are in one space dimension, the condition V • A = 
simply states that A is a constant. Moreover we have J(p) = — (l/2)\I/ / (p)p' + ty(p)E, 
where hereafter the apices denotes differentiation w.r.t. the macroscopic variable u. Note 
that condition ()6.3|) holds so that $ = U. Changing variables in (j4.11|) by introducing 
<p{u) = V(p(u)), u E [0,1]; we get 

U*(X) = ]- sup / du { - <p(u)ip'{u) 2 + AV(u) - V(") + 2XE^(u) 
2 ip Jo ^ 

where the supremum is over all positive <p such that <p(0) = ip , <p(i) = <fii and ip solves 

2^(u)V»'(m))' = ( - <p'(u) + 2[E + %(«))' 
with boundary conditions ip(0) = ip(l) = 0. The solution is 

^(u) = — log ^ + (E + X)u + A [ dv 



2 <Po Jo <p(v) 

where 

After elementary manipulations, the variational problem for U*(X) becomes 

4 y Jo 

The associated extremality condition, which determines the optimal profile, is 



2ip"{u)tp{u) - {<p'(u)) 2 - AE 2 ^{u) 2 = -a( - \ log ^ + E + a)* / dv -± 

\ 2 ipo ) IJ tp{v) 

with the boundary condition <p(0) = <p , <p(l) = <pi. 



i 1-2 

(6.13) 
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In the case E = it is not difficult to check that the solution of ()6.13j) is 



^ u ) = c{u + ^— x ) (u 



e x \ / y? e A 



1 — e -V V v 9 oe A — <£>i 
where C — — (l — e _A ) (</3o eA — V?i)- We then get 

fT(A) = - ^(0)] = i(l - e~ A ) (y, e A - </?i) 

In the case E ^ 0, the solution of f)6.13j) is instead given by 

tp(u) = c(e 2Eu -a) (e- 2Eu -b) 

where 

(p e 2E+x - (fx . 1 - e~ A_2jB 

a ■ 



ip e x — ip i 1 — e~ x 

and 

(1 -e- A )(^ e A -^i) 
(e 2E - 1)(1 -e" 2S ) 

Notice that this solution converges, as — > 0, to the solution with no external field. 
Plugging this solution into the variational formula for U*, we get that 

u*(\) = -\w{i)-v'm + \fa-<M = ^{r^b^^-^+^lh^- 1 )}- 

(6.14) 

We conclude this section showing that we may invert the order of limits in (J4.9|) for 
zero range models. For < x, y < N + 1, \x — y\ = 1, recall that we denote by J\f^' v 
the total number of jumps from x to y in the time interval [0, t}. For < x < N, let 
_ jyx,x+i _ ,x ^ e ^ e ^.Q^gj curren t over the bond (x, x + 1). Note that we are 
including the boundary bonds. 

Consider the limit as microscopic time t goes to infinity of the Laplace transform of 
the total current: 



A' 



e N (X) = -^m-logE^JexpjAiV-^Q^ 1 }" 

x=Q 

and notice that /n given by (|5.11jl is related to the Legendre transform of by 



f* N (\) = e^A^-log^M))}) • 



Notice furthermore that this expression does not depend on the initial condition rj N by 
ergodicity. 

Since two currents Q^ ,x+1 Qf v+1 differ only by surface terms, in the asymptotic t j oo, 
we may replace all currents by Q®' 1 and obtain that 



eiv(A) = ^-limJlogE^e^' 1 ] 

/V t^oo t 



N t^oo t 
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To compute the previous limit, we represent the zero range process in terms of inter- 
acting random walks. Let N be the total number of particles at time 0: iVo = ^2 x f] x {0)- 
We start labeling these particles. New particles entering the system at the boundary 
get new labels in an increasing order. Denote by X l (t) the position at time t of the 
2-th particle. X 1 performs a weakly asymmetric random walk on A^v with absorption at 
the boundary and mean g(l) exponential waiting times. X 2 does the same but its clock 
rates are affected by X 1 . If they occupy different sites, the X 2 -exponential has rate g(l), 
while if both occupy the same site, its exponential clock has rate g(2) — g(l) and so on. 
We need for this construction the function g to be increasing. Moreover the condition 
g(k + 1) — g(k) > a > guarantees that these random walks will hit the boundary with 
probability one. 

Let Wi (resp. Ui) the indicator function of the event that the i-th particle created at 
the left (resp. right) boundary is absorbed at N + 1 (resp. 0). Denote by N±(t) a Poisson 
process of rate tfo/2, ipi/2 which represents the entrance of particles at either boundary. 
With this notation, up to negligible terms in the limit t f oo, 

AT_(i) N+(t) 

Qt' 1 = Wi ~ Ui ■ 

i=l i=l 

Since the random variables Wi, Uj are independent, elementary computation shows that 
e^v(A) = ^{(<p /2)p N {e x - 1} + ( Vl /2)q N {e~ x - 1}} , 

where p^ = P[wi = 1] (resp. 1 — q^) is the probability that a random walk, absorbed in 
and N + 1, with transition probability p(x, x + 1) = e 2E ^ N / (e 2E ^ N + 1) = 1 — p(x, x — 1) 
starting from 1 (resp. N) is absorbed in iV + 1 (resp. 0). These probabilities can be 
explicitly computed. 

As N | oo, we get that 

i™^) = £ { r ^ E ( eA - 1 ) + ? #h(^- 1 )} 

which agrees with (j6.14|) . 

A. Supplement to Section OH 

We present here a derivation, at the heuristic level, of the law of large numbers and the 
large deviations, as — > oo, for the empirical density and the empirical current. Recall 
the notation introduced in Section El We have seen there that to prove the law of large 
numbers for the empirical measure and the current, we need to express the limit current 
J in terms of the density p. 

In the context of stochastic lattice gases this is done by assuming a local equilibrium 
state. Roughly speaking, this means that in a large microscopic region A around u, still 
infinitesimal macroscopically, the system has relaxed to the Gibbs state (with Hamiltonian 
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H) conditioned to YlxeA^x — \A\ir N (t,u). This assumption, which can be rigorously 
justified, [19], allows us to express the empirical current in terms of the empirical density. 
We next show how this can be done for the so-called gradient models. 

By standard computations in the theory of Markov processes we have that, [30, Lemma 
II. 2. 3], for a bond {x, x + e,}, 

Q x ^{t) = (l/2)iV 2 [ ds^+^-c^M] + M x ' x+e *(t), 

Jo 

where M x ' x+ej (t) are martingales with bracket 

{MX , X+ e %M y,y + e 3){t) = {1 /2 ) N 2 6 x ,ykj /^ [c x , x+ei (r/(s)) + C x+ei , x ( V (s))] ■ 



Let G be a smooth vector field as in ()3.2j) vanishing on dA. By definition of the martingales 

((J N , G)) T = ~ [ T dtJ2J2 G ^ X / N M C ^+eMt))-C, + eUvm , 

(A.l) 

where Ai^(G) is a martingale term. An easy computation, based on the explicit formula 
for the quadratic variations of the martingales M x,x+ei {t), shows that Ai^(G) vanishes 
as N — > oo. We next use definition (J2.2J) and Taylor expansion to write 

= 'I,.,,!'/) [1 + ^(x/N)] - c x+eu M [1 - ±Ei(x/N)] + 0(1/A 2 ) (A _ 2) 
= Kx + eM ~ + ^ K* +ei (v) + c x+ei> M]Ei(x/N) + 0(1/A 2 ) 

The gradient condition, see [30, II.2.4], holds if there exist local functions h^\rj), i = 
1, . . . , d, depending on the configuration rj around 0, so that for any i — 1, • ■ ■ , d 

c° x , x+ei (v) ~ c° x+ei ,M = h x l) (v) - h^M (A.3) 

where h$ is the function h$ evaluated on the configuration r\ translated by x. Let us plug 
the right hand side of (jA.2|) into (|A.1|) . By the gradient condition (jA.3|) we can perform 
a summation by parts on the first term. Note that there are no boundary terms since we 
assumed G to vanish on the boundary. We get, with a negligible error as N — > oo, 

{{jN > g))t w \w f dt E E ^(i) 

i=1 x (A.4) 
+ G,(t, x/N) [cl x+ei (v(t)) + c x+ei> Mt))]Ei(x/N)\ 
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Recall the definition of the functions S l \ x introduced in (|3.5|) . By the local equi- 
librium assumption mentioned above and the equivalence of ensembles from (|A.4|) we get 



((J n ,G))t « Yl f dt fdu^G^^d^i^^u^+G^^x^i^^u^E^u)} 
i=i Jo J A 

(A.5) 

Taking the limit iV — > oo, the empirical density ir N (t, u) converges to p(t, u), whereas the 
empirical current J^ N (t, u) converges to a vector field J(t, u). Equation (jA.5|) then implies 

J(p) = - l -D{p)Vp + X {p)E (A.6) 
where D and x are d x d diagonal matrices with entries Da(p) = -^d^(p) and Xu(p) = 

x (i) (p)- 

We now turn to a heuristic derivation of the large deviations principle ()3.8)) - (J3.9j) for 
the current. Recall the statement and the notation introduced in Section El 

In order to make the trajectory j typical, we introduce an extra weak time dependent 
external field F = (Fi, . . . , F d ) by perturbing the rates as in Section |2l namely 



Cxx+e 



(V) = c x , x+ M e N ~ lF ^ , c F x+ei M = c x+ei ,M e~ N - lp ^ . (A.7) 



TV F 

We denote by P the probability distribution of the perturbed process. Since these 
rates c F are the same as the rates of the original process with E replaced by E + F (cf. 
()2.2|) ). we have the following law of large numbers: 

where 

3 = J(p) + X(P)F = - l -D(p)Wp + X (P)(E + F) (A.8) 

and p satisfies the continuity equation (J3.10jh 

We now read this equation in the opposite direction: given the trajectory j we first 
solve ()3.10|) to get p, then we determine the external field F which makes j the typical 
behavior, namely 

F = x{pr 1 {j + \D{p)Vp)-E (A.9) 
By writing the original process in terms of the perturbed one we have 

d¥ N 

^ N [j N (t,u)^j(t,u), (t,u) E [0,71 x A ) =P^ F (— 

The large deviation principle ()3.8|) - (j3.9j) will follow, recalling (|A.9J) . once we compute 
the Radon-Nikodym derivative and show that on the event {J N ~ j} we have, with a 
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negligible error if N — > oo, 

dP N dF N,F 1 r T 

^ = -^~- Ni 2h it{FMF) (A ' 10) 

This equation can be interpreted, in analogy to the classical Ohm's law, as the total work 
done in the time interval [0, T] by the external field F . 

We shall need some basic tool from the general theory of jump Markov processes that 
we briefly recall, see e.g. [19, Appendix Al] or [4, Appendix A]. Let O be a countable 
set and consider a continuous time jump Markov process X t on the state space Q with 
generator given by 

Lf(v) = £ Kri)pM) [f(r)') - /(i?)] (A.ll) 
v'en 

where the rate A is a positive function on Q and p(rj, rj') is a transition probability. We 
consider also another process Xf of the same type with time dependent rate A F (?7, t) and 
transition probability p F (r],r]';t). Then, denoting by F Vo and P,^ o the distribution of the 
two processes with initial condition r] we have 



(X t ,te[o,T]) 

f ™ \ F (X T . , Ti)p F (x r , X T . ■, tA r T r F/ , A (A - 12) 

= ex P <^ v log v y- 1 ' , v ; - / dt X F (X t , t) - A I t > 

Itr \{x n _ 1 ) V {x n _ 1 ,x Ti ) j L 1 ; 1 UJ j 

where ro = 0, Xq = rjo, t*, % — 1, . . . , n is the time in which the process jumped from X Ti _ x 
to X Ti , and n is total number of jumps in the time interval [0, T\. 
For the process rj(t) with generator ()2.3|) we have for x, y in Z d , 

A^p^,^) = (1/2)^^(77) X(rj) = (1/2)N 2 J2 C *M 

For the process r](t) with rates (jA.7|) we have 

A^t^fo^+^t) = (l/2)N 2 c x , x+ei (v)e N ' lF ^ x/N) 
and a similar formula for \ F (r),t) p F (rj,a x+ei ' x ri,t) so that 



i=l x 

From ()A.12|) and the explicit expressions for the rates, we get that 

3 — 1 *C ^"x , x-\-e,j -\-ej : x 

^, x+e Mm(e N - lFit ' x/N) ~ l)+c x+ej ,Mt))(e~ N ~ lF ^ - l)} , 
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where t X)V are the jump times from x to y. Expanding the exponentials and recalling the 
definition of the empirical current, we may rewrite the previous expression as 



N d ((J N ,F)) T - - / dt^^Fiix/Nrffa^Wty-c^Mt))} 

x 1=1 

~ \ I dt ED^W')) + c ^Mt))} Fi{x/N,tf + 0(1/N) 

ZJo x i=X 



(A.13) 



where we let c X)V = if x, y £" A N . 

For gradient models condition (jA.3|) holds so that we can perform a summation by 
parts in the second term. Recalling the definition of the diffusion matrix D, the mobility 
X and the local equilibrium assumption, we can express the second term of the right hand 
side of (|A.13|) in terms of the empirical density. Since we are assuming J N j, we get 
that 

d^Y , ( r T i r T , 

log^r *N d {((j,F)) T + / dt(F,(l/2)D(p)Vp- X (p)E)-- dt(F, x (p)F)j . 

(A.14) 

which, by the choice of F in (|A.9[) . concludes the derivation of (|A.10|) . 

The rigorous derivation of action functional X requires some difficult estimates. In 
fact, while in the proof of the hydrodynamic limit it is enough to show that the local 
equilibrium assumption holds with a negligible error as — > oo, in the proof of the large 
deviations we need such an error to be o(e~ CN ). This can be proven by the so called 
super exponential estimate, see [19,21], which is the key point in the rigorous approach. 

Recall the dynamical large deviations principle for the density stated in (|1.4j) . The 
rate functional T is given by 

^[o,T](p) = \J^dt{VH{t),x{p{t))VH{t)) (A.15) 

where, given the fluctuation p, the external potential H = H(t,u) is chosen so that it 
vanishes at the boundary and 

dtp = V-(±D(p)Vp-x{p)[E + VH]) (A.16) 

which is a Poisson equation for H. 

We conclude this Appendix showing how the above result follows directly from the 
large deviation principle for the current. We fix a path p = p(t,u), (t,u) £ [0, T] x A. 
There are many possible trajectories j = j(t,u), differing by divergence free vector fields, 
such that the continuity equation (J3.1()|) is satisfied. The functional T[q,t\{p) can be 
obtained by minimizing X[ 0j t](j) among all such paths j 

F[o,T](p) = inf 1[o,T](j) ( A -!7) 

3 ■ 

Vj=-dtp 
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To derive the functional (jA.15|) we show that the infimum above is obtained when the 
external perturbation F introduced in (jA.9|) is a gradient vector field whose potential H 
solves (|A.16Jl . Let H be the solution of (JA.16J) and F as in (JA.9J) . we write 

F = VH + F (A.18) 

By the definition of H we get 

(VH, x(p)F) = -(h, V ■ j + V ■ (\d( p )V P - X {p)E - x(p)Vtf) ) = 

Hence 

%t] 0) = \l dt { ( V #> X(P)VH) + (F , x(p)F) } 

Therefore the infimum in (|A.17|) is obtained when F — 0, so that the functional defined 
in |A~T7jl coincides with (fA~T5jl . 
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